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SECTION  1 
INTRODUCTION 

The  understanding  and  analysis  of  three-dimensional  fluid 
behavior  around  a  surface  mounted  blunt  obstacle  (turret)  is  of 
considerable  importance  for  flight  applications.  The  prime 
motivation  to  solve  such  problems  has  been  to  investigate  the 
basic  phenomenon  of  separated  flow  and  to  determine  the  asr ->ciated 
forces  on  the  object.  In  addition,  the  dynamic  features  c-  the 
unsteady  flow  field  around  the  obstacle  are  of  interest  tc  :e 
designer . 

The  time  variation  of  the  separation  line  due  to  the  .<•  .ding 
of  vortices  for  the  flow  past  a  turret  produces  a  complex  and 
intriguing  flow  structure.  Also  the  appropriate  turbulence  model, 
ranges  of  Mach  number  and  Reynolds  number,  and  protuberance 
height  are  some  of  the  elements  in  parametric  space  which  restrict 
the  unified  approach  and  generality  of  the  solution.  The  experi¬ 
mental  and  computational  approaches  are  the  common  aerodynamics 
tools  to  investigate  complex  problems.  However,  during  recent 
years,  the  computational  approach  has  been  preferred  because  of 
its  lower  cost,  availability  of  faster  computers,  and  detailed 
description  of  the  fluid  physics.  In  this  frame  of  reference, 
the  present  work  addresses  the  "computational  flow  visualization" 
of  the  three-dimensional  separated  flow.  The  primary  objective 
is  to  provide  an  efficient  numerical  procedure  for  the  solution 
of  the  governing  equations  and  to  compare  the  computed  results 
with  the  previously  obtained  experimental  data. 

A  survey  of  the  literature  reveals  a  substantial  but  not  an 
exhaustive  effort  in  this  area.  The  experimental  results ^  ^ 
obtained  primarily  by  conventional  flow  visualization  techniques 
and  the  data  correlated  indicate  some  important  facts  about  the 
primary  separation  features  for  flow  past  small  protuberances. 

For  different  protuberance  heights,  the  available  flow  field  data 
is  so  sparse  that  its  interpretation  leading  to  a  general  conclusion 
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is  not  obvious.  To  acquire  the  finer  details,  it  is  felt  that  the 
three-dimensional  separated  flow  past  a  turret  requires  numerical 
analysis . 


Influenced  by  the  results  of  experimental  studies  and  a  search 
for  proper  criterion  for  three-dimensional  flow  separation  led  to 
a  simple  convincing  topological  approach1  .  This  approach  is 
based  on  the  hypothesis  that  vector  fields  of  the  skin  friction 
lines  and  external  streamlines  remain  continuous.  Taking  bodies 
of  revolution  at  various  angles  of  incidence,  for  example,  it  was 
possible  to  synthesize  the  flow  mechanism  in  a  rational  manner. 

We  will  be  using  some  of  these  topological  ideas  to  explain  the 
fluid  behavior  around  the  turret. 

Turning  to  our  present  effort,  for  subsonic  high  Reynolds 
number  flow,  past  a  hemispherically  capped  cylindrical  obstacle 
(the  turret),  the  horseshoe  vortex,  turbulent  boundary  layer 
characteristics  and  near  wake  structure,  need  careful  attention, 
particularly  from  a  numerical  simulation  point  of  view.  A  subset 
of  this  immense  problem,  unsteady  flow  oscillations  past  a  circular 
cylinder,  was  recently  studied  by  Shang^10^.  For  this  numerical 
experiment,  the  mesh  point  distribution  and  appropriate  boundary 
conditions  were  optimized  to  yield  the  time  mean  and  fluctuating 
properties  of  the  fluid.  This  two-dimensional  computation  provided 
encouragement  to  continue  further  and  to  analyze  the  large-scale 
organized  fluid  motion  around  a  three-dimensional  turret  at  a 
Mach  number  of  0.55  and  a  Reynolds  number  of  10.3  *  106/m.  At 
these  conditions  the  turret  height  is  much  larger  than  that  of 
the  turbulent  boundary  layer  on  the  surface.  The  aim  is  to  evaluate 
the  density  data  field  and  its  fluctuations  around  the  turret  and 
compare  the  results  with  the  experimental  data  of  deJonckheere, 
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SECTION  2 
ANALYSIS 

For  the  problem  under  consideration,  the  solution  of  unsteady 
compressible  Navier-Stokes  equations  becomes  necessary.  These 
equations  offer,  on  a  fine  scale,  a  good  insight  into  the  viscous- 
inviscid  interaction  and  the  separated  flow  structure.  For  the 
present  Reynolds  number  range,  the  flow  along  the  flat  plate  as  well 
as  around  the  turret  will  be  turbulent,  and,  thus,  adequate 
inclusion  of  turbulence  terms  is  needed.  In  the  absence  of  body 
forces  and  electromagnetic  effects,  the  governing  equations  can 
be  written  as 

+  v  •  (pu)  =  0 

a-(3£^)  +  V  •  (puu  -  t)  =0 

+  V  •  (q  -  u  .  V)  =  0 

where  the  heat  flux  q  is  given  by  the  Fourier's  law 
q  =  -K  V  T 

and  the  shear  stress  tensor  t  is  defined  as 

t  ^  j  =  -(u+e)(Def  ^  (u  +  e)V* 

In  order  to  optimize  the  numerical  resolution,  using  body  oriented 
coordinate  transformation,  these  equations,  in  mass  averaged 
variables,  can  be  written  in  the  following  convenient  form: 
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X 

p  +  2/3(p+e)(|H  +  +  |J) 

(6-1) 

xy 

=  xyx  =  ~(v+e)(|^  +  |5) 

(6-2) 

xz 

=  xzx  =  -<V+e>( H  +  |5) 

(6-3) 

y  “ 

p  +  2/3(p+e)(^+  |5  +  |D 

-  2(p+e).|5 

(6-4) 

yz 

=  xzy  =  “(p+e)(|5  +  If) 

(6-5) 

'z 

P+2/3<,,  +  .>(g+g+!5) 

-  2(p+e).|5 

(6-6) 

i 


(7-1) 


e  =  e.  +  (u2  +  v2  +  w2)/2.0  (8) 

This  system  of  equations  is  closed  by  a  perfect  gas  law, 
Sutherland's  viscosity  formula,  and  appropriate  eddy  viscosity 
model  as  a  supplementary  relation. 

The  fluid  motion  defined  by  the  aforementioned  set  of  equations 
is  to  be  analyzed  around  a  turret  (Figure  1)  mounted  on  a  flat 
surface.  Its  base  is  a  circular  cylinder  with  diameter  0.44  meter 
and  the  height  to  diameter  ratio  1.071.  One  of  the  sides  has  a 
cut-out  (window)  of  area  0.249  meter  *  0.201  meter  at  a  normal 
distance  0.194  meter  from  the  turret  axis.  The  top  of  the  turret 
is  a  hemispherical  dome  of  radius  0.223  meter  whose  center  is  on 
the  axis.  The  £ ,  n ,  correspond  to  normal,  axial,  and  tangential 
transformed  coordinates  of  the  turret,  respectively,  with  its 
origin  at  the  turret  base.  For  this  configuration,  the  experi¬ 
mental  results ^  ^  ^  were  available  for  various  locations  of  the 
turret  window  (look  angle) .  In  order  to  demonstrate,  a  typical 
flow  field  for  a  90  degree  look  angle  was  simulated. 

One  of  the  major  tasks  in  numerical  simulation  is  the  proper 
selection  of  mesh  system  in  the  computational  domain.  A  judicious 
choice  of  grid  points  helps  capture  the  important  flow  field 
characteristics  as  well  as  reducing  the  data  processing  time. 

After  much  study  the  following  choice  was  made.  The  turret  surface 
is  defined  by  62  x  30  points.  On  the  cylindrical  portion,  15 
points  were  used  in  the  axial  direction  stretched  exponentially 
from  the  base  and  15  points  were  uniformly  placed  on  the  hemi¬ 
spherical  dome.  In  the  tangential  direction,  62  equidistant 
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points  were  used  with  the  first  two  points  overlapping  to  permit 
the  use  of  fourth  order  pressure  damping  (Figure  2) .  Using  these 
surface  nodes  as  the  reference  points,  the  normal  coordinates  are 
then  described  by  30  exponentially  stretched  points  from  the  turret 
surface  extending  outward  to  ten  times  the  turret  radius  (Figures 
3  and  4).  The  finer  mesh  near  the  turret  and  flat  plate  will  help 
resolve  the  viscous  effects  while  the  coarser  mesh  helps  reduce 
the  computational  time.  Care  has  been  taken  that  no  mesh  point 
lies  on  the  turret  axis;  however,  a  sufficient  number  of  mesh 
points  is  provided  to  resolve  the  flow  field.  For  this  mesh 
arrangement,  the  derivatives  of  coordinate  transformation  can  be 
evaluated  through  the  Jacobian,  J. 


This  Jacobian  of  the  coordinate  transformation  is  a  non-singular 
determinant  and  allows  a  one-to-one  correspondence  from  the 
physical  to  the  transformed  space. 

In  any  numerical  algorithm,  the  boundary  conditions  play  a 
crucial  role  in  both  the  accuracy  and  stability  of  the  numerical 
scheme.  The  initial  and  boundary  conditions  associated  with 
this  system  of  equations  are  straightforward  and  simple.  Initially, 
the  turret  is  assumed  to  be  immersed  completely  in  a  uniform 
stream.  The  upstream  boundary  conditions  are  merely  the  free 
stream  value,  whereas  at  the  downstream  side,  the  gradient  of 
properties  is  assumed  to  vanish.  At  the  solid  contours,  the  no 
slip  conditions  for  velocity,  ortho-isobaric  condition  for  density 
and  isothermal  conditions  are  imposed;  while  near  the  turret  axis, 
an  interpolation  is  adopted  to  eliminate  calculation  on  this 
singularity  (Figure  5) .  Thus,  for  the  independent  variables, 

0  <  (,  n»  c  <  1,  in  the  computational  domain,  the  conditions  are 
set  in  the  following  forms. 
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Figure  2.  Mesh  Point  Distribution  on  the  Turret. 
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Initial  Conditions 


U(0,e,n,C)  =  Uro(c,n,c) , 

Boundary  Conditions 

Upstream:  U(t,£,l,e)  =  U(t,£fi,sQ), 


'l  <  50  *  c2 


Downstream:  3  „ _  ,  N  _  n 

U  (t,  c  #  n  f  Cq)  =  0 


At  the  Solid  Contours,  n-0  u  =  v  =  w  =  0 

T  =  T  =  1.0584T 
w  <*> 


P,.  =  P„ 

W  00 


w 


At  the  Axis  U(t,C,l,c)  =  F (U (t , 5 , 1 , cQ) , Cj. , C2 ' AL) 


(9) 

(10) 

(ID 

(12-1) 

(12-2) 

(12-3) 

(13) 


where  ^  and  are  the  extremum  of  the  left  half  of  the  computa¬ 
tional  domain  exposed  to  the  free  stream  and  AL  is  the  arc  length 
subtended  by  two  consecutive  points  along  the  turret  axis  to  the 

center  of  the  hemisphere.  The  significant  variation  in  the  turret 

[12] 

surface  curvature  and  strong  vortical  interactions  require 

careful  selection  of  the  turbulence  model.  Once  the  turbulence 
model  is  identified,  the  system  of  equations  is  closed  and  the 
computation  can  be  initiated. 

Eddy  Viscosity  Model 

The  closure  of  the  system  of  equations  is  achieved  by  using 
a  Cebeci-Smith  type  eddy  viscosity  model  and  by  assigning  a 
turbulent  Prandtl  number  of  0.90.  The  specific  turbulence  model 

[13] 

used  is  a  simplified  version  due  to  Baldwin  and  Lomax 
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In  the  inner  region: 

=  p  •  (K^  *L*D^)  2  •  |  to  |  (14) 

where  to  is  the  vorticity  of  the  flow  field.  In  the  present 
formulation,  the  Van  Driest 's  damping  factor  is  given  as 


D 


1 


1 


and  the  scaling  length  L 
normal  from  the  surface. 

In  the  wake  region. 


Vp*/|fc>w[ 


exp  (--\  /  —  -  -  -  •  L/26)  (15) 

is  the  distance  measured  along  the  outward 


we  use  the  modified  Clauser's  law. 


eo  =  0-00168'p»*Umax*6i  (16) 

where  6*  is  the  kinematic  displacement  thickness  of  the  wake  in 
the  plane  of  the  turret  shoulder.  This  thickness  is  evaluated  at 
a  downstream  distance  about  six  times  the  turret  radius. 


r*  _  f  u2+v2+w2\ 

11 V 


dc 


The  effective  eddy  viscosity  is,  thus,  taken  to  be 


e  =  min (e  .  e . )  . 

o  1 


(17) 
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SECTION  3 

NUMERICAL  PROCEDURE 

The  present  investigation  utilizes  the  well  established  vector 
solver  '  which  adopts,  basically,  the  MacCormack's  explicit 
finite  difference  scheme .  Owing  to  the  large  memory  availability 
and  high  data  processing  rate  (the  CPU  time  per  grid  point  per 
time  iteration) ,  the  CRAY-1  computer  was  choosen  for  our  computation. 
The  intricate  interplay  between  its  capacity  and  numerical  scheme 
provided  an  efficient  and  powerful  tool  to  analyze  this  three- 
dimensional  problem  with  acceptable  engineering  accuracy. 

We  used  the  solver  to  initialize  the  planes  or  "pages"  in  the 
(n-C)  plane  and  to  sweep,  for  each  time  step,  in  the  5 -direction, 
thus  minimizing  the  data  flow  from  and  to  the  disk  or  mass  storage 
devices.  The  vector  length  for  this  case  was  62,  which  is  also 
the  longest  index  array  of  the  three-dimensional  mesh  system. 

The  computational  facilities  were  utilized  with  remote  access  via 
long  distance  dial-up  and  the  SAMNET  network  link.  The  data 
processing  rate  of  6.8  *  10-5  seconds  was  achieved  for  the  present 
analysis.  Most  of  the  available  memory  was  required  to  perform 
the  time  dependent  calculations. 

The  finer  mesh  sizes  in  the  normal  direction  of  the  turret 
and  flat  surface  were  0.357  x  10  2meter  and  0.345  x  10  2meter, 
respectively,  and  the  angular  displacement  between  nodes  was  six 
degrees  (6°) .  These  mesh  sizes  seem  to  provide  satisfactory  time 
step  increments  during  the  computation.  Further  refinement  could 
result  in  the  stiffened  behavior  of  the  transformation  metric 
element  near  the  dome  and/or  much  smaller  time  step.  A  conservative 
choice  of  the  Courant-Fredrich-Lewy  (CFL)  condition  was  implemented 
to  have  optimum  time  step  with  stable  numerical  evolution. 

Once  the  initial  phase  of  computation  was  completed  and  some 
anticipated  periodicity  in  the  flow  characteristics  was  observed, 
the  evaluation  of  the  root  mean  square  (RMS)  density  fluctuation 
was  attempted.  For  turbulent  compressible  flow,  separating  the 
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time  mean  and  fluctuations  of  density,  the  squared  density 
fluctuation  over  a  characteristic  time  period  was  calculated  as 


<p 


->  2 


> 


(18) 


where  is  the  principal  time  period,  equal  to  0.0114  seconds 
and  is  a  function  of  Strouhal  number,  turret  diameter  and  free 
stream  velocity.  The  computation  of  RMS  density  for  all  grid 
points  demanded  three  times  the  additional  available  computer 
memory  which  was  not  feasible  at  this  stage.  Therefore,  this 

particular  analysis  was  confined  to  a  part  of  the  computational 

.  y  Y 

domain  from  a  turret  height  -  =  0.4  to  r  =  1.7.  A  total  of 

K  K 

9880.0  seconds  CPU  time  of  CRAY-1  computer  were  required  (using 
a  CRAY  FORTRAN  compiler  (CFT)  version  1.09)  to  calculate  the  entire 
flow  and  the  RMS  density  fluctuations  over  one  principal  period. 
Almost  "all"  the  memory  resources  of  the  computer  were  required 
to  process  the  three-dimensional  separated  flow. 
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SECTION  4 

RESULTS  AND  DISCUSSIONS 

This  section  contains  two  main  areas  of  discussion.  The  first 
of  these  will  deal  with  the  comparison  of  the  computed  and  experi¬ 
mental  results ^ ^ ^ .  This  permits  the  validation  of  our  computa¬ 
tional  procedure  as  well  as  providing  guidelines  for  the  subsequent 
analysis.  The  rest  of  the  discussion  is  concentrated  on  probing 
into  the  different  aspects  of  three-dimensional  separated  flow 
structure. 

The  pressure  coefficient  distribution  at  the  turret  shoulder 
and  45°  meridian  plane  angle  (MPA)  are  presented  in  Figure  6. 

The  variable  THETA  is  the  azimuth  angle  measured  clockwise  from 
the  oncoming  free  stream  direction.  Thus,  the  cavity  side  of  the 
turret  is  located  in  the  range  0°  to  180°  and  the  smooth  side  is 
between  180°  and  360°.  We  will  adhere  to  the  same  definition 
throughout  our  discussion.  The  comparison  shows  good  agreement 
with  the  experimental  results.  The  different  curves  represent  the 
variation  in  pressure  coefficient  during  one  principal  time  period. 
At  the  turret  shoulder,  where  the  cylindrical  and  hemispherical 
surfaces  meet,  the  computed  mean  CP  (Figure  6a)  is  0.9918  for 
the  forward  stagnation  point  (AZA  =  0°),  whereas  the  measured 
value  is  0.91.  The  front  edge  of  the  turret  window  experiences 
a  sudden  expansion  of  the  flow  and  a  small  separation  zone  resulting 
in  a  steep  drop  in  pressure  (AZA  =  63°)  and  a  rapid  pressure 
recovery  within  20°  in  the  downstream  direction.  The  smooth  side 
does  not  witness  such  abrupt  changes  and  gives  way  to  minimum 
pressure  value  at  AZA  =  279°.  For  the  45°  MPA,  on  the  hemispherical 
surface  (Figure  6b) ,  the  CP  distribution  exhibits  relatively  smooth 
variation.  It  also  hints  that  the  presence  of  the  turret  window 
has  very  small  effects  for  this  meridian  plane. 

This  argument  is  further  confirmed  by  Figure  7,  which  shows 
the  instantaneous  pressure  distribution  contours  of  different 
levels  around  the  turret  surface.  A  zone  of  very  low  pressure  can 
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Figure  6.  Comparison  of  Pressure  Distribution. 
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Figure 


be  observed  at  the  turret  window  spread  over  10°  azimuth  angles 
extending  from  slightly  above  the  flat  surface  on  to  the  hemisphere 
up  to  37°  MPA.  The  low  valued  pressure  contours  are  mostly  con¬ 
fined  to  the  rear  side  of  the  turret  for  the  azimuth  angles  80° 
to  300°.  At  the  turret  top,  the  averaged  computed  value  is 
-1.27  which  is  lower  than  the  experimental  value  -1.0.  It  may 
be  recalled  that  the  turret  top  is  in  the  vicinity  of  60  clustered 
mesh  points,  none  of  which  is  at  the  top.  Our  pressure  coefficient 
value  is  the  mean  of  the  values  at  all  these  points.  This  dif¬ 
ference  may  be  attributed  to  the  grid  point  arrangement  and  the 
improvement  in  this  direction  is  one  of  our  future  goals. 

A  significant  insight  into  the  separated  flow  is  offered  by 
the  limiting  streamline  patterns,  which  are,  in  the  limiting  case, 
the  skin  friction  lines  on  the  flat  surface  around  the  turret. 

In  Figure  8,  we  present  a  sequential  development  of  such  lines 
that  agree,  at  a  particular  instant  of  time,  qualitatively  to  the 
work  of  Sedney  and  Kitchens^.  The  primary  separation  line 
located  at  twice  the  turret  radius  upstream  of  the  forward  stag¬ 
nation  point  and  the  attachment  lines  are  quite  clear  and  distinct. 
However,  the  secondary  separation,  appearing  on  and  off  near  the 
primary  separation  line  during  the  process  of  computation  are 
difficult  to  detect.  These  lines  merge  into  the  primary  separation 
line,  very  close  to  the  plane  of  symmetry  indicating  that  it  lifts 
off  the  surface  very  quickly.  This  observation  will  be  further 
confirmed  when  we  discuss  the  flow  field  in  the  cross-sectional 
plane.  The  primary  separation  line  that  encircles  the  turret 
front  and  extends  downstream  in  the  asymmetric  plane  finally 
becomes  a  leg  of  the  horseshoe  vortex. 

The  asymmetry  in  the  plane  induced  by  the  turret  window  con¬ 
tains  some  spectacular  features  in  the  downstream  direction.  The 
vortical  motion  in  the  near  wake  is  produced  by  the  adverse  pres¬ 
sure  gradient  at  this  separation  point  and  rear  turret  window 
edge.  These  separation  points  on  the  surface  together  with  rear 
stagnation  point  qualify  for  the  half  saddle  point  (Hunt,  et  al.^). 
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The  turret  window  introduces  an  additional  node  and  saddle  point 
close  to  the  surface.  Accounting  for  the  node  of  separation  in 
the  downstream  direction,  saddle  point  at  the  plane  of  symmetry 
at  the  primary  separation,  and  regrouping  them,  the  following 
topological  rule  for  streamlines  on  a  two-dimensional  plane  cutting 
the  three-dimensional  body  is  satisfied 


<en  +  1  V>  -  lIS  +  1  V>  -  -1  <19' 

where  is  the  sum  of  the  nodes  and  is  the  sume  of  the  saddle 
points.  This  confirmation  provides  a  sound  reasoning  for  coupling 
the  analytical,  experimental,  and  computational  studies. 

The  asymmetry  in  the  flow  is  further  exemplified  in  Figure  9, 
showing  the  local  Mach  number  distribution  for  the  azimuth  angle 
of  90°  (cavity  side)  and  270°  (smooth  side) .  Close  to  the  flat 
plate  (Figure  9a) ,  the  cavity  side  records  lower  Mach  number. 

The  recovery  is  only  77%  of  the  free  stream  value  at  the  outer 
boundary  which  is  less  than  5  percent  of  its  value  at  the  smooth 
side  outer  boundary.  This  may  be  associated  with  the  strong 
vortical  flow  near  the  rear  edge  of  turret  window  and  the  dominant 
viscous  effects.  An  altogether  different  situation  is  observed  at 
the  turret  shoulder  (Figure  9b) ,  where  the  Mach  number  attains  its 
peak  values  of  0.83  and  0.62,  very  close  to  the  body  and  asymptoti¬ 
cally  approaches  the  free  stream  value  within  1.5  turret  radii. 

The  experimental  data^^'^®^  also  endorses  these  values. 

In  Figure  10  we  present  the  RMS  density  fluctuations,  normalized 
by  free  stream  density  for  the  azimuth  angle  90°.  Each  curve  repre¬ 
sents  the  density  fluctuations,  at  a  fixed  turret  height,  along  a 
normal  direction.  The  peak  value,  0.0405,  happens  to  be  at  the 
turret  window  base.  As  we  move  upwards,  this  peak  value  shows 
a  decreasing  trend  and  then  oscillates  between  0.009  and  0.003 
(Figure  11).  The  corresponding  peak  value ^ ^  (experimental)  is 
0.034,  19  percent  lower  than  our  computed  result.  This  experi¬ 
mental  value  was  deduced  from  the  hot  wire  data  using  an  empirical 
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Figure  9.  Mach  Number  Distribution. 
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Figure  10.  RMS  Density  Fluctuations  for  the  90°  Azimuth  Angle. 


relation  and  ignoring  the  higher  order  pressure  fluctuations. 

The  far  field  fluctuations,  which  are  small,  might  not  be  of  major 
concern,  but  in  the  near  field  the  turret  base  is  most  sensitive 
to  such  fluctuations.  This  fact  is  further  elaborated  in  Figure 
12,  giving  RMS  density  fluctuation  contours  at  the  base,  middle, 
and  top  of  the  turret  window.  It  is  interesting  to  observe  that 
these  fluctuations  are  grouped  at  two  different  locations  and  the 
peak  value,  0.073,  occurs  at  AZA  =  246°.  These  fluctuations  of 
high  intensity  slip  off  towards  the  downstream  direction  and  cluster 
between  AZA  =  100°  and  210°  with  a  peak  value  of  0.04  at  AZA  = 

190°.  The  density  in  the  upstream  direction  remains  almost 
invariant  and  does  not  give  rise  to  any  substantial  fluctuations. 

Figure  13,  14,  15  give  a  glimpse  of  the  flow  field  structure 
in  the  cross-sectional  plane  (azimuth  angle  0°  -  180°)  of  the 
computational  domain.  It  presents  the  instantaneous  velocity  field, 
Mach  number,  and  density  contours  in  a  sequential  manner.  The 
significant  upstream  influence  is  confined  close  to  the  flat  plate 
up  to  2.5  times  the  turret  radius.  The  velocity  field  is  sug¬ 
gestive  of  an  interesting  feature:  along  the  front  stagnation 
line  on  the  cylindrical  surface,  the  flow  is  almost  divided  into 

Y 

two  parts.  At  —  =  0.3543,  the  contribution  of  the  axial  velocity 
component  is  neglibible  compared  to  the  tangential  and  normal 
components.  The  accelerating  flow  at  the  hemispherical  dome,  with 
a  maximum  Mach  number  of  0.87  at  12°  MPA  upstream,  is  separated 
at  47°  MPA  in  the  downstream  direction.  The  reverse  flow  in  the 
near  wake  region  and  at  the  turret  shoulder  and  their  time  depen¬ 
dent  variation  provide  the  mechanism  for  vortex  shedding  and  the 
relative  movement  of  reattachment  point. 

Complementary  to  this  cross-sectional  plane,  the  flow  field 
distribution  in  the  plane  parallel  to  the  flat  plate  offers  a 
qualitative  comparison  to  flow  past  a  cylinder.  In  this  plane, 
up  to  the  turret  shoulder,  the  fluid  phenomenon  shows  a  similar 
trend.  As  a  typical  example,  the  velocity  field  distribution, 

Mach  number,  and  density  contours  are  presented  at  two  different 
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Figure  15.  (Concluded)  Density  Contours  in  the  Cross 
sectional  Plane  at  Different  Times. 
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locations  (Figures  16  through  18,  at  —  =  0.3543  and  Figures  19 

K 

through  21,  at  turret  shoulder).  The  turret  window  contributes 
an  additional  amplification  to  the  process  of  vortex  shedding  in 
the  downstream  direction.  At  the  turret  shoulder,  the  maximum 
local  Mach  number  attained  was  0.93  at  AZA  =  63°  and  0.90  at  AZA  = 
270°.  These  values  are  much  lower  than  that  for  the  corresponding 
two-dimensional  case^^  ,  where  transition  and  even  shock  wave 
formation  is  reported.  Also,  the  separation  point  has  advanced 
further  downstream  towards  the  rear  stagnation  point  at  AZA  =  99° 
and  243°,  respectively.  Even  though  the  separation  point  location 
is  a  function  of  time,  it  is  observed  that  for  a  particular 
instant  of  time,  the  locus  of  the  separation  point,  the  separation 
line,  changes  its  course  considerably  along  the  turret  height. 

This  illuminating  fact  is  further  revealed  in  Figure  22. 


This  Figure  shows  an  instantaneous  limiting  streamline  pattern 
on  the  turret  surface.  Near  the  stagnation  point,  a  saddle  point 
type  singularity,  a  foci,  node  of  attachment  along  the  turret 
height  can  be  easily  identified.  The  turret  window  includes  an 
additional  small  separation  zone  and,  hence,  introduces  a  node 
and  saddle  point.  These  singular  points,  together  with  the  nodes 
at  the  forward  and  rearward  stagnation  points,  give  another  con¬ 
firmation  to  the  topological  criteria.  For  the  skin  friction  lines 
on  a  three-dimensional  body,  B,  connected  simply  to  a  plane  wall. 


P,  the  sum  of  nodes,  en 
the  following  relation 


and  the  sum  of  saddle  points, 

[6] 


Eg  satisfy 


<en  - 


ZSJP+B 


The  asymmetric  separation  lines  start  just  above  the  foci,  at 
the  turret  base,  climb  up  towards  the  hemisphere  and  then  meet  at 
the  rear  stagnation  line  at  30°  MPA.  The  separation  effects  on 
the  turret  top  are  mostly  confined  to  63°  MPA.  On  the  cylindrical 
surface,  near  the  turret  shoulder,  the  separation  points  show  a 
tendency  to  move  upstream  indicating  that  strong  three-dimensional 
flow  exists. 
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20.  (Concluded)  Mach  Number  Contours  Around  the 
Turret  Shoulder  at  Different  Times. 


Finally,  these  limiting  streamline  patterns  offer  a  precise 

description  of  the  downstream  flow.  The  situation  is  similar  to 

flow  past  a  hemisphere  cylinder  at  90°  incidence.  The  separation 

line  around  the  turret  passes  through  the  singular  points, 

indicating  a  close-type  of  separation.  This  registers  an  inacces- 

r  19 1 

ability  of  the  outer  flow  to  the  near  wake  region  .  This 
interpretation  follows  logically  from  the  steady  state  separated 
flow.  However,  the  periodic  mass  exchange  due  to  unsteady  effects 
is  evident 


50 


SECTION  5 
CONCLUSIONS 

An  efficient  numerical  procedure  is  successfully  accomplished 
for  the  three-dimensional  separated  flow  around  a  turret.  An 
adequate  description  of  flow  structure  is  obtained  and  an  endorse¬ 
ment  of  a  topological  criteria  is  presented.  For  a  free  stream 
Mach  number  of  0.55  and  a  Reynolds  number  of  10.3  *  10 6 /meter, 
good  comparison  is  established  between  our  computed  results  and 
the  experimental  data  for  both  the  pressure  distribution  and  levels 
of  RMS  density  fluctuations. 
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